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We investigate the dynamics of two interacting diffusing particles in an infinite effectively one 
dimensional system; the particles interact through a step-like potential of width b and height 4>o an d 
are allowed to pass one another. By solving the corresponding 2+1-variate Fokker-Planck equation 
an exact result for the two particle conditional probability density function (PDF) is obtained for 
arbitrary initial particle positions. From the two-particle PDF we obtain the overtake probability, 
i.e. the probability that the two particles has exchanged positions at time t compared to the initial 
configuration. In addition, we calculate the trapping probability, i.e. the probability that the two 
particles are trapped close to each other (within the barrier width b) at time t, which is mainly of 
interest for an attractive potential, 4>o < 0. We also investigate the tagged particle PDF, relevant 
for describing the dynamics of one particle which is fluorescently labeled. Our analytic results are 
in excellent agreement with the results of stochastic simulations, which are performed using the 
Gillespie algorithm. 

I. INTRODUCTION 

As recent advances in manufacturing methods drives device sizes toward the nanorange, the understanding of how 
interactions between diffusing entities affect dynamics is becoming increasingly important^. Situations where diffusing 
molecules interact strongly are also of importance in biological systems^. 

The interaction between diffusing particles can be of either attractive or repulsive nature. For repulsive interactions, 
a particularly prominent example is that of single-file diffusion, i.e. the diffusion of identical particles which interact 
via a hardcore repulsion (the interaction potential energy is plus infinity, so that the particles cannot pass each other) 
in one dimension. For single-filing systems the particle order is thus conserved over time (t) resulting in interesting 
dynamical behavior for a tagged particle. For instance, in contrast to ordinary diffusion for which the mean square 
displacement ((x-r — xt,o) 2 } is proportional to t, for single file diffusion the mean square displacement of a tagged 
particle is proportional square root of time, {{xq- — xj-.o) 2 ) oc t 1 / 2 for long times in an infinite system with a fix 
particle concentration^^; the probability density function (PDF) the single-file of the tagged particle position is 
G aussian2&£&2 i 1 1 1 1 1 1 2 . For attractive interactions an especially well-studied example is that of reaction-diffusion 
system^, where often the particles are assumed to annihilate each other upon encounter, i.e. the potential energy 
between particles is assumed to be minus infinity. 

Although much work has been dedicated to interacting diffusing particles interacting via infinite (negative or 
positive) potentials, to our knowledge, much fewer studies consider finite potentials. In Ref. [13 the problem of 
diffusion of N particles on two coupled linear chains was studied. Similarly, in Refs. [H and[l6| diffusion of spherical 
particles in a cylindrical geometry, where the cylinder radius was large enough to allow passage of particles, were 
studied; in Ref. [TH molecular dynamics simulations were done using a Lennard Jones interaction between particles 
and in Ref. [TH a Monte Carlo simulation using a hard sphere interaction was performed. Recently, the dynamics of a 
tagged particle in a system consisting of particles interacting through screened repulsive Coulomb interactions in one 
dimension was investigated^, however only through stochastic simulations. In this study we derive analytic results for 
diffusing particles interacting through finite potentials: we solve analytically the problem of diffusion of two particles 
interacting via a finite-sized potential of finite height in one dimension, for arbitrary initial particle positions. Our 
results generalize the single barrier results of Ref. [l8| (who solved a similar problem by Laplace transform techniques) 
to arbitrary initial particle positions. 

This paper has the following organization: In Sec. [Til we state the problem under consideration and formulate the 
relevant equations. In Sec. IIIII we provide the solution of the equations for the two particle conditional probability 
density function (PDF). In Sec. HVI we use the two particle PDF to obtain the overtake probability, i.e. the probability 
that the two particles at time t have exchanged positions compared to the initial configuration. In Sec. [V]we calculate 
the trapping probability i.e. the probability that the two particles are trapped close to each other (within the barrier 
width) at time t. In Sec. IVII we obtain the PDF for one of the particles being tagged. We compare our analytic 
results to stochastic simulations using the Gillespie algorithm and find excellent agreement. Finally, in Sec I VIII we 
give a summary and outlook. 
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FIG. 1: Cartoon of the problem considered in this study: two point particles are diffusing in an infinite one dimension system. 
Particle j has coordinate Xj, and initial coordinate Xj t o (j — 1,2). The particles interact via a step-potential of width b and 
height <f>o, see Fig. [2] 
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FIG. 2: Phase space region for the system, the darker area, |xi — xs\ < 6/2, corresponds to a region where the potential is of 
height 0o, the potential width is b. In the lighter shaded areas the potential is zero. The system is assumed to be infinite, i.e. 

— CO < Xl,X2 < oo. 



II. PROBLEM DEFINITION 



We consider a system with two interacting point particles diffusing in an infinite one dimensional system. The 
point particle problem considered here can be transformed into a problem of finite-sized interacting particles using a 
similar mapping as given in Ref. \vX A cartoon of the system we have in mind is depicted in Fig. [T] The particles 
have coordinates x = (£1,2:2) and initial positions xq — (a?i,o> ^2,0)- We assume that the particles interact through a 
step- like potential of width b and height <f>a, i.e. the potential is 

<&(x\,x 2 ) = for xi — x 2 < —b/2 (region 12) 

= 4>q for \x\ — xa| < b/2 (region 1 — 2) 

= for xi - x 2 > b/2 (region 21). (1) 

The phase space is depicted in Fig. [2j where the darker area corresponds to the region of non-zero potential. For 
0o > we have a barrier, whereas for 0o < the potential is of a short-range attractive nature. 

The spatial distribution of the particles as a function of time is contained in the two particle conditional PDF 
V(x, t\xo), which gives the probability of particle 1 being in the interval [x\, x\ + dxi] and particle 2 in [x 2 , x 2 + dx 2 ] 
at time t given that they initially (at time t = 0) were at positions x\q and £2,0 respectively. Since inside each of 
the three regions the potential energy landscape is flat (the force is zero) this quantity is governed by the 2 + 1 (two 
coordinates and time) variable Fokker-Planck (Smoluchowski) equation 
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where 7 = 12, 1 — 2 and 21 indicates phase-space regions, and D is the diffusion constant for the particles. At the 
boundaries between regions we have the following condition s 19 ' 20 

flux J continuous at \xi — x 2 \ = b/2 

cxp(— f3^(x))Vj(x, t\xo) continuous at \xi — x 2 \ = b/2, (3) 

where (3 = 1/(&bT), and ks is the Boltzmann constant and T the temperature. We assume that the particles start 
in region 12 initially, i.e particle 1 starts to the left of particle 2. The initial condition then becomes: 

Vl 2 (x,t — > 0\x ) = S(X! -Xlfl)6(x 2 -x 2 ,o), 

Vi- 2 (x,t ^ Q\x Q ) = V 21 (x,t^0\x )=0, (4) 
where 6(z) denotes the Dirac delta-function. 



III. TWO-PARTICLE PROBABILITY DENSITY FUNCTION 

In order to solve the equations specified in the previous section we make a variable transformation to the ccnter- 
of-mass position X and relative coordinate s according to: 

X = \ ( x i + x 2 ), 
s = x\ — x 2 . (5) 

Eqs. |[2J) and (j4|) then becomes (leaving the argument corresponding to the initial positions implicit) 

&P^(X,s,t) / x 5 2 , ns d- 



V 12 (X,s,t^0) = S(X-X )S(s-s ), 
Vx- 2 {X,s,t^Q) = V 2X {X,s,t ^0) = 0, (6) 

where Xq = (xi.q + x 2> q)/2 and sq = xi.q — x 2 ,o ( s o < ~b/2) and the effective diffusion constants 

x D 
D x = — , 
2 ' 

D s = 2D. (7) 

The equations above express the fact that the relative coordinate s diffuses with a diffusion constant 2D, whereas the 
center-of-mass coordinate X diffuses with a diffusion constant D/2. The boundary conditions, Eqs. give rise to 
the four equations: 

dT 12 (X,s,t) i d7V 2 (AVs,t) , 

Wa ls =- b/2 ~ ds U =- fa/2 ' 

&P 1 - 2 (X,s,t) i dV 2 i(X,s,t) t 

7\ s=b/2 — 7{ S=6/2) 

os 1 os 1 

which expresses the continuity of flux and 

V 12 (X, s = -b/2, t) = cxp(-^-)7V 2 (X, a = -b/2, t), 

exp(^=)Pi- 2 (^,s = 6/2,*) = V 21 {X,s = b/2, t), (9) 
k B T 

which corresponds to a detailed balance-like condition 2 -!, at the boundaries. The boundary conditions above are 
different from the corresponding problem in quantum mechanics^ 2 - where the wavefunction and the derivative of the 
wavefunction are continuous - the origin of this difference is due to the fact that in the quantum mechanical problem 
the potential enters in the equation of motion (the Schrodinger equation), whereas in the Fokker-Planck equation it 
is the force which enters. 



(8) 



4 



Eqs. ([6]), (8J and ([9]) allow a product solution of the form 

V 7 (X,s,t) = V x (X,t)P°(s,t), (10) 
where (the boundary conditions involve only the relative coordinate s) 

For the solution for ^(s, i) we take a general function that satisfies the diffusion equation in each of the regions and 
that satisfies the boundary conditions: 

V s ( Sl t) = r ^. e -m)t e -iQso x [ C + e iQs + c-e-*^] for s < -6/2 

./-oo 27r 

x [C+_ 2 e lQs + Ci_ 2 e~ iQs ] for |s| < 6/2 

x C&e" 2 * for s > 6/2. (12) 

The prefactors are dependent on Q in general. Inserting Eq. (fT2|) into the equation of motion [see Eq. gives 
the dispersion relation: 

E(Q)=D S Q 2 . (13) 

We proceed by setting C-jt> = 1 (we will show that with this choice the initial condition is satisfied). Inserting Eq. 
(fT2"]) into the boundary conditions Eqs. (jSJ) and @ produces 4 equations for the four unknowns C^ 2 , C^_ 2 and Cj/_ 2 
and C21. Solving this set of equations gives: 
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cr_ 2 (Q) = -ni-^) 



1 _ p2 e 2iQb ' 
1 



^ _ p2 e 2iQb 



CUQ) = ^ + m-F) 1 _p 2e2iQb , (14) 



where we introduced an effective "reflection" coefficient 



F 



_ = tanh[^]. (15) 



Notice that for the case of a barrier- like potential, 4>o > 0, we have < F < 1, where F = corresponds to the 
absence of the barrier and F — > 1 corresponds to an infinite barrier. For the case of an attractive potential, 4>o < 0; 
we have —1 < F < 0, where F — > — 1 corresponds to an infinite potential well. Combining Eq. (I12[) . (Ti"3"l) an( i (U3J) wc 
have: 

n 2 (M) = *(s - s , *)+ F{T(-[s + s + 6], t)-T(-[s + s -&],*)}, 
Pi-a(M) = [l-F)?(s-s ,t)-F(l-F)?(-[s + s -b],t), 
V s 21 {s,t) = (l-F 2 )T( s - So ,t), (16) 



where 



and 



*( a , t) = / ^^_ ( D-D 3 Q 2 t iaQ = exp (_^L) (17) 



/J riU 6 "^- (18) 



5 



The above expression for T(a,t) can be explicitly evaluated: using 1/(1 — z) = X^^Lo z ™ (valid for \z\ < 1) we have 



t(m = r^f> 2 - 



o D s Q't e iaQ(a+2bn) 



E 3 ^ 1 , (a + 26n) 2 

= / (S^TJ^— 4D^ } - (19) 

i.e. T(a,t) is a sum of shifted Gaussians weighted by powers of the reflection coefficient F [sec Eq. (p~5|) ] . The full 
solution to the problem is specified by Eqs. pH|) (fTf]) , P^|) . (JTTJ) and (TH))) , where X and s are related to x\ and x 2 
using Eq. ([5]). Finally, returning to our original coordinates we find that the two-particle conditional PDF becomes: 

Vi2(x,t\x ) = Q(xi - xi,o)Q(x 2 - x 2 , ) 

+j p ^ F 2n {g(x 1 - £2,0 + g - 6n)0(a!2 - a?i,o - ^ + 5n ) 

n=0 

-Q(xi - x 2 ,o - - - bn)Q(x 2 - Xi t o + - + bn)} 

OO 

V X -2{S,t\S ) = {l-F)J2 F2n {Q(xi-xi,o + bn)Q(x 2 -x 2 ,o-bn) 



n=0 



-FQ(x% - x 2 ,o ~\~ bn)Q(x 2 - £i,o + ^ + bn)} 



T 21 (x,t\x ) = (l-F 2 )Y,F 2n Q(x 1 -x 1 , + bn)Q(x 2 -X2,o-bn), (20) 



n=0 



where 



^-(C^-^S'- (21) 



It is a straightforward matter to show, using 8(x — a) = limt_,o cxp[— (x — a) 2 / (4Dt)]/y4vrDi, that the result above 
satisfies the initial condition specified in Eq. © [and that, therefore, indeed taking Cf 2 = 1 was the correct choice]. 
In the absence of a barrier, F — ► 0, we find V(xi, x 2 , t) = Q(xi — Xi$)Q{x 2 — x 2 .o) as it should. For the case of an 
infinite barrier, F — ► 1, we have: V\ 2 {x\, x 2 , t) = Q(a;i — Xi,o)<9(£2 — £2,0) + Q(a:i — x 2 ,o + b/2)Q(x 2 — 0=1,0 — 
and Vi- 2 (xi, x 2l t) = 7 7 2i(^i, ^2, t) = in agreement with the result for hardcore interacting particles of linear size 
b/2*£ For the case of an infinite potential well F — > — 1, the solution for region 12 takes the form: Pi 2 (xi,x 2 ,t) = 
Q(x\ — xi_o)Q(x 2 — x 2 fl) — Q(xi — x 2 .o + b/2)Q(x 2 — x± t o — b/2); this result agrees with previous results for two vicious 
walkers (the random walkers kill each other upon encounter)] 23 ' 24 

IV. OVERTAKE PROBABILITY 

In this section we calculate the overtake probability W(t), i.e. the probability that particle 1 is to the right of 
particle 2 at time t. We have: 



W(t) = J ds:i / dx 2 \ rcgion 2iV(x u x 2 ,t). (22) 
Changing coordinates to X and s, see Eq. (JSJ, we get 

/oo />oo 
dXV x (X,t) / dsT s 21 {s : t) 
-oo Jb/2 



= I(l-F 2 )E^erfc[^±i§^], (23) 



V4DH 

n—0 v 



where erfc(z) = 1 — erf(z) is the complementary error function, erf(z) = (2/ 'y/n) Jq dy exp(— y 2 ) is the error function 2 ^, 
and so = %i,o ~ £2,0 as before. We have above used the fact that J^ oo dXP x (X,t) = 1 [see Eqs. (fTTj) ] together with 
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FIG. 3: Overtake probability W(t) as a function of time t for different barrier heights <f>o- The solid curves are the analytic 
result, Eq. ([23}, and the dashed curves are the results of stochastic simulations. The following parameters were used xi t o = —3b, 
£2,0 = 3b, where b is the barrier width. The simulation was performed using the Gillespie algorithm on a lattice with M — 5000 
lattice sites, averaged over n = 10000 ensembles, and with a lattice spacing a = 0.0256. 

Eqs. ([TP]). |H]), (TTT]) an d (QT)]). The result given in Eq. agrees with the result derived in Ref. QJ (using 

a Laplace-space formalism), where the problem of passage of one particle (since the center-of-mass coordinate is 
integrated out above, W(t) is effectively a one-dimensional quan tity) across one and two barriers in finite and infinite 
one dimensional systems were considered - however in Ref. [l8| the initial position were taken to be right at the left 
border (i.e. so = —b/2), the result given in Eq. (|23p thus generalizes the one barrier (infinite system) result in Ref. 
[H to general initial condition. We point out that ( 0o <-» — </>o is equivalent to F *-> —F) since W(t) contain only 
even powers of F the overtake probability is invariant under the reversal of sign of the potential <fio <-> — </>o; it thus 
takes equally long times to pass a barrier of height <fio as it takes to pass a trap of height — 4>q. Using the fact that 
crfc(0) = 1 wc find that for long times t — > oo we have W(t) — > 1/2 as it should. 

In Fig. [3] we illustrate the overtake probability as given in Eq. (|23[) . We notice that an increased potential height 
leads to decreased probability for overtaking. Also, due to the time it takes for the particles to approach each other 
through diffusion, there is an initial time r c before the probability of overtaking becomes appreciable (even for zero 
potential); a simple estimate gives t c ~ (^2,o — xi o ~ b) 2 /(2D s ). In Fig. [3] we also compare the analytic results to 
that of a stochastic simulation using the Gillespie algorithm (see Appendix [A]) , and find excellent agreement. 



In this section we calculate the trapping probability T(t), i.e. the probability that the two particles are at a distance 
closer than b to each other at time t. We have 



V. 



TRAPPING PROBABILITY 




(24) 
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FIG. 4: Trapping probability T(t) as a function of time t for different potential heights (j>o. The solid curves are the analytic 
result, Eq. (|25[l . and the dashed curves are the results of a stochastic simulation. Except for the differently chosen values for 
0o we used the same parameters as in Fig. [3] 



Again, changing coordinates to X and s, see Eq. ([5]), we have 



/oo 
dXV x (X 7 t) / dsV?_ 2 (s,t), 
-oo J-b/2 



(l_Z)g F 2« {(1 + F)erf[ 



71=0 



b(2n + 1/2) - sp 



CTf[ b(2n-l/2)- S0 ] _ Fai[ b(2n + 3/2)- S0 ]] _ 



\/4D s t 



s/4DH 



(25) 



For long times t — > oo we have that the trapping probability approaches zero T(f ) — > 0; this is due to ergodicity [each 
point in our system is assigned a probability proportional to exp(— f3<f>o)] combined with the fact that our potential 
well is of finite width (= b) and connected to an infinite system. In Fig. [¥]we illustrate the trapping probability given 
by Eq. ([25|k we find excellent agreement with stochastic simulations. We notice that a more negative potential depth 
leads to increased probability for trapping as it should. Similarly to the overtake probability, due to the time it takes 
for the particles to approach each other through diffusion, there is an initial time r c before any significant number of 
trapping events has occurred. 

A possible experiment testing the predictions in this section would involve, for instance, two fluorescent molecules; 
the two molecules interact through a potential and assuming that when the molecules are within some distance b from 
each other the total fluorescence get quenched or enhanced one could directly detect a binding event between the 
two particles as an increase or decrease in total fluorescence (fluorescence measurement are here assumed ensemble 
averaged over thermal noise, with fixed initial particle positions). In such an experiment it may be more convenient 
to, rather than obtain T(t) from experimental data, use time-integrated fluorescences, i.e., measure the time-averaged 
trapping probability as given by % v (t) = (1/t) f T(t')dt'. With the help of Eq. (JUJ) we straightforwardly obtain: 



f1 — F~) °° 

%At) = ^r^E i?2n {( l + F ) G ( 



2 



n=0 

b(2n - 1/2) - s , 



&(2re+l/2) - s 



V4DH 



FC M2n + 3/2)-s 



(26) 



where we defined a function G(z) = 2zexp(— z 2 )/y / 7r — 2_z 2 crfc(z) + crf(z). In Fig. \5\ we illustrate the time-averaged 
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trapping probability given by Eq. (|26p for different potential heights as a function of t. Experimental measurements 
of T(t) or 7^ v (£) would provide detailed information about the nature of the interaction potential between the two 
particles, i.e. of <f>o and 6, by comparison to our analytic expressions. 

Finally as a simple check of the results in Sees. llVl andlV l we calculate the first passage time density for two vicious 
walkers (F = —1). The survival probability S(t) (i.e. the probability that particle 1 and 2 has not met at time t) is 
given by S(t) = 1-W(t)-T(t), which for F = -1 becomes S(t) = erf[-(s + b/2)/y/4D s t], using Eqs. {33) and (35). 
The first passage time density is then -dS(t)/dt = [\s + b/2\/\/4nD s }t- 3 ^ 2 cxp[-(s + b/2) 2 /{AD a t)\, which agrees 
with the standard result^l for the first passage time problem for a diffusing particle (diffusion constant D s ) starting 
a distance sq + b/2 from a perfectly absorbing wall, as it should. 

VI. TAGGED PARTICLE PROBABILITY DENSITY 

By integrating out x 2 we can obtain the tagged particle PDF pi(xi,t\xo) for particle 1. We have: 

/oo 
dx 2 V(x 1 ,x 2 ,t) 
-OO 

r'Xx+b/2 



/ dx 2 Vi2{xi,x 2 ,t) + / dx 2 Vi~ 2 (x l7 x 2 ,t) 

Jxx+b/2 Jxi-b/2 

/Xi—b/2 
dx 2 V 21 (x 1 ,x 2} t), (27) 
-OO 
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and, using Eq. (|20p . we find the explicit expression 

Pl(xi,t\x ) = Q(xi - Xifi)R{xx - X 2 fi + -) 

°° 6 
+F^F 2n [Q{x 1 - x 2 .o + x - bn)R{xi - x lfi + bn) 

n=0 

—FQ(xi - x 2 ,o — - - bn)R(x 1 - x xfi + b + bn)] 

oo , 

+(1 -F)J2 F 2n {^Q{ Xl - x lfi + bn)R{ Xl - x 2 ,o + ■=- hn ) 

—FQ(x\ — x-ifi — — — bn)R(xi — xi + bn) 
+Q(?i - xi fi + bn)[(l +F) — FR( Xl - x 2fi bn)}}, (28) 

where we introduced the function: R(a) = crfc[a/(V ADt)]/2. The corresponding result for the PDF for particle 2 is 
obtained by the replacements x\ —X2 and x\fi <-> — X2.0 in Eq. (|28|l . In the absence of a barrier, F — > 0, we find 
/Oi(o;i, t\xo) = Q{x\ — 2:1,0), i-c the tagged particle PDF is that of an independent diffusing particle, as it should. 
For the case of an infinite barrier, F — > 1, we have: pi(x\, t|xo) = Q(xi — Xi t o)R(xi — X2.0 + 6/2) + Q(x\ — X2fi + 
b/2)R(xi — xi,o) in agreement with the result for hardcore interacting particles of linear size &/2 i 12 i 26 The tagged 
particle PDF calculated here can experimentally be investigated, for instance, by fluoresccntly labeling of one of the 
particles. 

VII. SUMMARY AND OUTLOOK 

We have in this study solved exactly the problem of diffusion of two particles interacting via a step-like potential of 
height 0o and finite width b in an infinite one-dimensional system. In particular, from our exact analytic expression 
for the two particle probability density function, we obtained the overtake probability (the probability that the two 
particles has exchanged positions at time t), Eq. the trapping probability (i.e. the probability that the particles 

are at distances closer than 6/2 at time t), Eq. (|25[) . and the tagged particle probability density function, Eq. ([251) . 

For the case of a positive potential, <fto > our results arc of interest for the diffusive dynamics of repulsive 
particles. For instance, our results (and future extension to many particles) will be of interest for understanding the 
recent simulation results dealing with the diffusive dynamics of a tagged particles in a system of charged particles (of 
the same charge) interacting through screened Coulomb interaction^. 

For the case of an attractive potential 4>q < 0, our results should be of interest for understanding reaction-diffusion 
systems, where the reaction potential between the interacting species is of finite height and width. 

It remains a future challenge to extend the results of this study to many particles; in particular, it will be interesting 
to see how the t 1 ^ 2 (see Introduction) scaling of the mean square displacement of a tagged particle in a system of 
hardcore interacting particles and how the dynamics of reaction-diffusion systems arc modified as the potential height 
</>o is made finite. 
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APPENDIX A: STOCHASTIC SIMULATIONS 

Stochastic simulations using the Gillespie algorith m 27 i 28 i 29 is a technique well suited for generating stochastic 
trajectories for interacting particles systems. Briefly, (similarly to Ref. [l2.24.30.3l]) we consider hopping of two 
particles on a one-dimensional lattice, with lattice constant a. The number of lattice sites is denoted by M and 
chosen sufficiently large so that the ends of the lattice are not reached. The dynamics is governed by the 'reaction' 
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probability density function 

P(t, /Lt) = fc„ ex P(" 2 k ^ ^ A1 ) 

where r is the waiting time between jumps and k^ are the corresponding jump rates. There are four jump rates for the 
two particle system considered in this study: the rate for particle 1 jumping to the left (right) is k\ (fe); similarly we 
denote by &3 (fej,) the left (right) jump rates for particle 2. For the case that the two particles are at a distance smaller 
or larger than half the barrier width 6/2, the two particles diffuse independently and we set k\ = k 2 = fc 3 = fc 4 = fcf ree - 
If two particles are separated by a distance 6/2 we set reduced jump rates according to k% = k% = fcf ree exp(— P<f>o). 
We generate a stochastic time series through the steps: (1) place the particles at their initial positions; (2) From 
the PDF given in Eq. (|A1[) we generate the random numbers r (waiting time) and \i (what particles to move and 
in what direction) using the direct method^ 7 -; (3) Update the positions Xi (i = 1,2) of the particles, the time t and 
the rates ki,k2,k% and k^ for the new configuration and return to (1). This procedure produces a stochastic time 
series Xi(t). Steps (l)-(3) are repeated n times (n ensembles) in order to obtain the ensemble averaged overtake and 
trapping probabilities or the tagged particles PDFs (which are the entities calculated in the main text). The ensemble 
averaged results of a Gillespie time series is equivalent to the solution of a master equation incorporating the rates 
given abov o 27 i 28 . In the limit a —> with fixed diffusion constant, D = fcf reo a 2 , the master equation approaches the 
diffusion equation as specified in section HT1 
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